Acceleration Statistics as Measures of Statistical Persistence of Streamlines in 

Isotropic Turbulence 
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We introduce the velocity V s of stagnation points as a means to characterise and measure statis- 
tical persistence of streamlines. Using theoretical arguments, Direct Numerical Simulations (DNS) 
and Kinematic Simulations (KS) of three-dimensional isotropic turbulence for different ratios of 
inner to outer length scales L/n of the self-similar range, we show that a frame exists where the 
average (V s ) = 0, that the r.m.s. values of acceleration, turbulent fluid velocity and V s are related by 
La'/u' 2 ~ (Vs/u'){L/rj) 2/3+q and that Vi/u' ~ {L/r/) q with q = -1/3 in Kolmogorov turbulence, 
q — — 1/6 in current DNS and q = in our KS. The statistical persistence hypothesis is closely 
related to the Tennekes sweeping hypothesis. 
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Accelerations are central to fluid flow. Acceleration 
statistics are central to turbulence dynamics, turbulent 
mixing and transport, and also important where the tur- 
bulence controls droplet growth, chemical reactions, com- 
bustion and other processes (e.g. UQl). A key statistic 
is the acceleration variance. Its value and scaling with 
Reynolds number are essential to stochastic Lagrangian 
models and to Lagrangian probability density function 
models of turbulent diffusion if these models are to in- 
corporate finite Reynolds number effects ■ In this letter 
we argue that the Reynolds number scaling of the turbu- 
lence acceleration statistics, in particular the acceleration 
variance, are also important as measures of the statistical 
persistence of streamlines. 

The persistence of streamlines is central to turbulent 
pair diffusion 0, 0, IS 01 It may a l so be important for 
one-particle two-time Lagrangian turbulent statistics . 
If streamlines are persistent enough in time, fluid ele- 
ment trajectories approximately follow them for signif- 
icantly long times. Hence, initially close fluid element 
pairs typically separate when they encounter a region of 
highly curved diverging streamlines around a straining 
stagnation point. However, streamlines and their persis- 
tence are not Galilean invariant. This issue motivated 
to suggest that the dependence of turbulent diffu- 
sion on streamline structure fully emerges only in the 
frame of reference where streamline persistence is max- 
imised in some statistical sense. The assumption that 
such a frame exists, and that the statistical persistence 
of the streamline topology in this frame is long enough 
to leave its defining imprint on turbulent diffusion, they 
termed 'statistical persistence hypothesis'. This is a pow- 
erful hypothesis as it leads to predictions supported by 
DNS and KS such as 7 = 2d/D s (see H @) which 
relate the Richardson pair-diffusion exponent 7 (a La- 



grangian quantity) in d-dimensional turbulence to the 
fractal dimension D s of the spatial distribution of stagna- 
tion points (an Eulcrian quantity) in the privileged frame 
of the statistical persistence hypothesis. For isotropic 
turbulence, [f| argued that this frame is the one where the 
mean fluid velocity vanishes. The purpose of this letter is 
to mathematically formulate the statistical persistence of 
streamlines as well as the statistical persistence hypoth- 
esis (statistical persistence of streamlines maximised in 
a well-chosen frame) and to argue that, in isotropic tur- 
bulence, they are both increasingly valid for increasing 
Reynolds number. 

It is the persistence of highly curved diverging stream- 
lines around straining stagnation points which is impor- 
tant for turbulent diffusion, and it is therefore critical to 
consider the speed with which stagnation points move in 
space. Given an arbitrary frame of reference, the fluid 
velocity it at a stagnation point s(t) at time t vanishes, 
i.e. u(s, t) — 0, and remains so for as long as this stag- 
nation point exists. Hence, during the stagnation point's 
life-time, = j^u{s, t) — 4^ + V s ■ Vit at position s and 
time t, and = 4f is the stagnation point velocity. The 



fluid acceleration is defined as a = 
positions x and times t. Setting x 



a 



-V s ■ Vm 



^ + m ■ Vm at all 
s, we obtain 
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Note that relations (|TJ and (0) hold in any frame of ref- 
erence. What changes with frame is the number and po- 
sitions of stagnation points (where these relations hold) . 
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These relations are kinematic and are the fundamen- 
tal link we use here between accelerations and the persis- 
tence of stagnation points measured by V s . We formulate 
the concept of statistical persistence of streamlines on the 
basis of QJ. Assuming, as is reasonable in fully developed 
isotropic turbulence that the kinetic energy dissipa- 
tion rate per unit mass, e, scales as u' 3 / L (where u' is the 
turbulence velocity fluctuations r.m.s., and L is an outer 
length-scale of the turbulence) and that the velocity gra- 
dients in scale with the inner length-scale ij and the 
small-scale velocity u„ ~ (erj) 1 ^ 3 , then the acceleration 
r.m.s. a' is related to the r.m.s. V s ' of V s by 

V. ~ <*% (3) 
where r n ~ i]/u rn or equivalently 

La'/u' 2 ~ M/u'XL/rj) 2 ' 3 . (4) 

Strictly, a' is the acceleration r.m.s. over all stagnation 
points, but it is also equal to the acceleration r.m.s. over 
the entire field because Galilean transformations leave a' 
unchanged even though they cause the r.m.s. statistics 
to be calculated over different ensembles of points. Our 
DNS and KS calculations of a' and V£ support this view. 
The Kolmogorov scaling of a' (where a' is determined by 
e and rf) is La'/u' 2 ~ {L/ij) 1 / 3 and has been corroborated 
(using rj ~ (z^/e) 1 / 4 where v is the kinematic viscosity) 
by laboratory measurements of acceleration statistics in 
isotropic turbulence [|| in the range 900 < Re\ < 2000, 
where Re\ is the Taylor- length-based Reynolds number. 
This scaling implies V//u' ~ (L/^) -1 ' 3 which means that 
the movement of stagnation points (characterised by V,') 
is, statistically, increasingly slower than that of fluid ele- 
ments (characterised by u') as the Reynolds number in- 
creases. In this statistical sense, streamline curvature 
around stagnation points becomes increasingly persistent 
and fluid element trajectories can follow the curvature 
of these streamlines for an increasingly significant time. 
This interpretation assumes, as we confirm in our DNS 
and KS studies of isotropic turbulence which we report 
below, that in the frame where the mean flow vanishes, 
the average of V s over the entire flow is 0, i.e. (V s ) = 0. 
In other inertial frames, (V s ) is not zero but proportional 
to the velocity U of the frame relative to the one where 
the mean flow vanishes, thus reducing the persistence of 
streamlines in such other frames. (Stagnation points in 
these frames F correspond to points with fluid velocity U 
in the frame Fq where (V s ) = 0. The average acceleration 
over these points in frame Fo is (a\U), i.e. the average of 
a conditional on points where the fluid velocity is U, and 
is equal to the average acceleration over the correspond- 
ing stagnation points in frame F which follows from (JTJ 
and is given by — (V s ■ Vu). The DNS study which we 
report below indicates that (a\U) oc U for all Reynolds 
numbers tested, thus suggesting (V s ) oc (7, a propor- 
tionality relation confirmed by our DNS with a negative 
proportionality coefficient.) In conclusion, a frame exists 
where the statistical persistence of streamlines is max- 
imised and where the statistical persistence hypothesis 



is valid provided that V' a is much smaller than u'\ this 
frame is the one where the average of V s is zero, hence 
the one where the mean flow vanishes as argued by |(| . 

DNS studies of three-dimensional isotropic turbulence 
support the non-Kolmogorov scaling La'/u' 2 ~ (L/rf) 1 / 2 
in the range 40 < Re\ < 230 (see ^(J an d note that 
e ~ u' 3 / L is approximately valid in such simulations with 
such Reynolds numbers y| where, also, r\ ~ {v 3 /t) 1 ^) 
and therefore imply V s /u' ~ (L/^) -1 / 6 in that range. 
This also means that stagnation points become increas- 
ingly persistent as Reynolds number increases, but at a 
rate slower than that of Kolmogorov turbulence. Defin- 
ing an exponent, q, by 

La'/u 12 ~ (L/ V ) 2 / 3 +« ~ Re\ +3q/2 (5) 

(where Re\ ~ (L/rf)3 has been used), it follows from Q 
that (0 is equivalent to 

Vl/u>~{L/ V y~Rel«' 2 ; (6) 

q > corresponds to absence and q < corresponds to 
presence of statistical persistence of streamlines (in the 
sense that frames exist where streamlines are statistically 
fairly persistent and a privileged frame also exists — the 
one where (Vs) =0 — where this statistical persistence 
is maximised and the statistical persistence hypothesis 
holds). In the remainder of this letter we use DNS and 
KS to confirm and further explore the behaviour of V s , 
the Reynolds number scalings of La'/u' 2 and V s /u' and 
the relation between these scalings via (0} . An advantage 
of formula © is that it can be used to calculate V s from 
Eulerian snapshots of DNS and KS velocity fields without 
having to track the motion of stagnation points. 

We use DNS data of non-decaying homogeneous 
isotropic three-dimensional incompressible turbulence 
generated by a standard spectral code (the magnitudes 
of a few low wavenumber modes are kept constant) 
with grid resolution of about 2rj (except for the high- 
est Reynolds number where it is about 377) and compute 
stagnation points and their velocities V s in instantaneous 
velocity fields for Rex ranging from 57 to 250 (see [H for 
the DNS numerical scheme and parameters). We use the 
Newton-Raphson method to find all stagnation points. 
This is an iterative method and requires starting points, 
which have been taken over the DNS field's entire do- 
main at points separated by the numerical grid width. 
For every Reynolds number, we calculate (V s ) in various 
frames F, (a\u) in frame Fo and also a', Vg, the integral 
length-scale L and 77 (we calculate — Vp + v\7 2 u to ob- 
tain a. If we wrongly approximate a = — Vp, then we 
find (a\u) = 0). Our results show for all our Reynolds 

numbers, that (a\u) ~ ~Jl( u ~ ( see E a ); 

note that we observe components (fli\Uj) to be only 
if i ^ j) and that (V s ) = only when U = (see 
Fig- Gib)). We also plot La'/u' 2 and V^/(a'T n ) as func- 
tions of L/r] (see Fig. Pc)). In agreement with previous 
DNS studies (see and references therein) we obtain 
La'/u' 2 ~ (L/ry) 1 / 2 (i.e. q = —1/6) and confirm J3J| over 



3 




Ui/v! Ui/u' L/rj 



FIG. 1: Statistics of acceleration and stagnation point velocity. DNS. Solid triangles, Re\ = 57; open squares, 83; solid squares, 
120; open circles, 180; solid circles, 250. (a) Conditional average of acceleration in the frame Fo. (b) Average stagnation point 
velocity in the moving frame with U relative to Fo- (c) R.m.s. acceleration and r.m.s. stagnation velocity; here we plot results 
of ten different snapshots from each simulation. L/rj fluctuates in time which explains the small scatter around its average. 



our DNS range of Reynolds numbers. By the way, from 
{a\U} = - (V s • Vit) and (a\U) « -j^U one might ex- 
pect (V s ) (x U with a constant of proportionality of de- 
creasing absolute value for increasing Re\, as we indeed 
observe. What is unexpected (but is left for future inves- 
tigation) is the rate of this decrease; Fig. Cfb) suggests 
(V a ) = -&(£/r?)- x / 3 Uwhere& is a positive dimensionless 
universal constant. At any rate, our DNS suggest that 
(V s ) /v! -» and V'Jv! as Re\ -> oo. From Q, 
this implies that as Re\ — ► oo, stagnation points tend to 
become non-moving zero-acceleration points. 

To gain further insight into (V s ) and the validity of re- 
lations @, JSJ and iJBJ and how they depend on specific 
properties of the underlying flow, we now study synthetic 
velocity fields, namely KS, where the spectrum and the 
time-dependence can be modified at will 0, 01 IS E3- 
Such a study cannot be carried out with DNS where 
the spectrum and time-dependence of the flow are de- 
termined by the Navier-Stokes dynamics and cannot be 
tampered with. An additional advantage of KS is that 
the Lagrangian statistics it produces compare well with 
various DNS and laboratory results (see [l^] and refer- 
ences therein) . Finally, because of the dramatic decima- 
tion in number of modes, it is possible with KS to explore 
scalings with L/rj up to extremely large L/rj values (here 
10 6 ) which are out of reach of current DNS. 

In our KS we use three-dimensional turbulent-like ve- 
locity fields of the form (see 0] f° r fuller details) 



vectors A n and B n are determined from the energy spec- 
trum E(k n ) prescribed to be of the form 



u 



N k 

£ 

n=l 



A n cos(fc„ • x + U! n t) + B n sin(fc„ ■ x + u n t) 



where Nk is the number of modes, the directions and ori- 
entations of A n and B n are chosen randomly and uncor- 
rected with the directions and orientations of all other 
wave modes but perpendicular to fc„. The distribution 
of wavenumbers is geometric, specifically fc„ = |fe n | = 
feil.07 n_1 . The velocity field is incompressible by con- 
struction, and also statistically stationary, homogeneous 
and isotropic as shown by 0, Q . The amplitudes of the 



E(k) 



3(p- l)u r 
2(L/2tt)p- 



(7) 



in the range 2ir/L — k\ < k < fcjy fc = ^ir/r], and E(k) = 
otherwise. The ratio L/rj is increased by increasing Nk- 
Following 0,0], we set w n = Xu'k n , for different values 
of the dimensionless parameter A. 

From Jy we can derive a generalised form of Q within 
the framework of KS by assuming that the velocity gra- 
dients in (JJJ scale with the inner length-scale rj and the 

p— i 

small-scale velocity u„ ~ u (r]/L)^. This leads to 
La'/u' 2 ~ (Vg /u')(L/r?i)^ E ~ which generalises (@J and to 
the statement that La'/u' 2 ~ (L/r))^ E+q is equivalent 
to Vg/u' ~ (L/rj) q which generalise © and JHJ. 

KS runs with p = 5/3, values of L/rj ranging between 
10 and 10 3 , and A = 0,0.5,5, all lead to (a\u) = and 
to (V s ) = in frame Fo. In frames F, (V s ) = — U. In 
KS, V s is uncorrelated with Vtt so that (V s ■ Vtt) = 
(V s ) ■ (Vu) = in agreement with = (a\U) = 
— (Vs • Vtt). Runs with p = 1.4, 5/3, 1.8, values of L/rj 
ranging between 10 and 10 6 , and A = 0, 0.5, 5 also lead to 
La'/u' 2 ~ {L/rj)^ (see Fig. Ufa)) and therefore q = 0, 
and to Vg/u' — c where c is a dimensionless constant 
independent of L/rj, which confirms that q = (see 
Fig- |3[b)). (In all our KS cases, 5 x 10 6 starting points for 
the Newton-Raphson method are chosen over the same 
volume L 3 where L s is the L corresponding to the largest 
L/rj tried, keeping rj constant.) 

The result (a\u) — reflects the lack of dynamics 
and related lack of correlations between Fourier modes in 
KS. Turbulence dynamics seem to generate restoring ac- 
celerations which are anti-correlated with velocities (see 
Fig. [Ha)). 

The constant c turns out to be an increasing function 
of A (from our simulations, c oc A) and, as expected, c = 
for A = 0. Hence, in our KS isotropic turbulence where 
q = 0, the statistical persistence of streamlines is mea- 
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FIG. 2: Statistics of acceleration and stagnation point velocity. KS. (a) R.m.s. acceleration scaling with L/rj. For p — 5/3; 
crosses, A = 0; circles, 0.5; triangles, 5.0 For A = 0.5; asterisks, p = 1.4; squares, 1.8. (Data sets for these p are shifted with 
respect to the vertical axis by factor 10 for p = 1.4 and by 0.1 for 1.8.) The slopes of the lines are (3 — p)/2. (b) Variation of 
r.m.s stagnation point velocity with L/rj. Triangles, A = 0.5; crosses, A = 5.0. 



sured by c and is a direct reflection of the unsteadiness 
parameter A. In spite of q being different from the Kol- 
mogorov value —1/3, Richardson exponents 7 = 2d/D s 
(in particular 7 = 3 for p = 5/3) are observed in KS but 
only for small values of A 0, Q thus confirming the view 
that these exponents require some statistical persistence 
of streamlines to be realised 0- 

Finally, it is worth recalling the Tennekes sweeping 
hypothesis 01 which states that the dissipative eddies 
are swept past an Eulerian observer in a time much 
shorter than the time scale characterizing their own dy- 
namics. The statistical persistence hypothesis, the valid- 
ity of which we confirm in this letter, states that there 
exists a frame where (V s ) = and V/. <C u'; from 
it therefore follows that a' <C u'/r n . In this sense, the 
accelerations are small, which is a way to restate the 
Tennekes sweeping hypothesis. Indeed, the time needed 
for dissipative eddies to sweep past an Eulerian observer 
is r]/u' which is therefore much smaller than u n /a', the 
time which characterizes the dynamics of these eddies. 
This is the Tennekes sweeping hypothesis derived from 
the statistical persistence of streamlines. Alternatively, 
stagnation points mark regions of the flow where there 
is no sweeping. According to 0], statistics which are 
taken so as to remove the sweeping effect depend only 
on the small-scale dynamics, and this must therefore be 



the case of V,'. In Kolmogorov turbulence, the scaling 
of these small-scale dynamics is determined by e and v 
which therefore implies V' s ~ (e^) 1 / 4 in agreement with 
q = —1/3 and the statistical persistence hypothesis. In 
the present KS, however, sweeping of smaller-scale tur- 
bulence by larger-scale eddies is absent, and stagnation 
points correspond to regions where zeros of Fourier modes 
congregate. These Fourier modes move together with ve- 
locity w{k)/k. Hence, stagnation points typically move 
with that same speed and Vg ~ ui(k)/k ~ \v! in agree- 
ment with our KS results q — and c cx A. 

The statistical persistence of streamlines seems to be a 
reformulation and generalisation of the Tennekes sweep- 
ing hypothesis in terms of streamline topology and its 
persistence. Also, in conjunction with the kinematic re- 
lation Kolmogorov dimensional analysis implies the 
existence of a coherent flow structure, namely the per- 
sistent multiple-scale stagnation point structure of the 
turbulence. Elsewhere |14| we argue that the mean life- 
time of stagnation points is of the order of the integral 
time-scale of the flow. 
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